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We examine the effect of First Forbidden (FF) decays on /3-decay neutrino spectra by performing 
microscopic nuclear structure calculations. By analyzing the FF decay branches of even-even nuclei 
we conclude that FF decays may be responsible for part of the missing neutrinos in the so called 
’’Reactor Neutrino Anomaly”. Further calculations and more experimental data are needed for a 
firm conclusion. 
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I. INTRODUCTION 

The ” reactor anti-neutrino anomaly” is the observation 
that the average of the experimentally determined reac¬ 
tor anti-neutrino flux at reactor-detector distances less 
than 100 m accounts for only 0.946 ± 0.023 of the theo¬ 
retical expectation [H-Q . One of the explanations for this 
anomaly is that standard neutrinos that carry the weak 
nuclear charge can oscillate into a “sterile” neutrino that 
does not contain a weak nuclear charge and escapes de¬ 
tection. This is one of the few types of experiments that 
could be sensitive to the sterile neutrinos. It is important 
to verify that the shape of the calculated anti-neutrino 
spectrum is correct. 

In the effects of various corrections to the expected 
neutrino spectra were examined, but it was found that 
these corrections could not explain the anomaly. How¬ 
ever, in 0 only allowed decays were analyzed and the 
First Forbidden(FF) decay was left out. In Q the con¬ 
tribution of FF decays were examined by assuming some 
virtual FF branches in the decay, and it was found that 
the results could match the magnitude of the missing 
flux. However, the actual FF decays are much more com¬ 
plicated than the form assumed in Q. In this work, we 
will exam the effect of realistic FF decay from micro¬ 
scopic nuclear structure calculations, and make an esti¬ 
mation of their effect on the neutrino spectra. 

For the neutron-induced fission accumulated yield dis¬ 
tributions one finds two peaks for the fission products. 
One of these is centered near ®‘*Sr and another is centered 
near ^'^°Xe. Both or these regions could have comparable 
amount of FF branches. For the region of ^^°Xe calcu¬ 
lations with the Quasi-Particle Random Phase (QRPA) 
method were carried out Q, and reasonable agreement 
between experiment and theory for half-lives and logit 
values were obtained. Also shell-model (SM) calculations 
for some nuclei in the Xe region can be carried out in a 
moderately large basis with previously derived Hamilto¬ 
nian. Thus, in this work we will focus on applying the 
QRPA and SM methods in this Xe region to investigate 
the effect of FF beta decay on the shape of the neutrino 
spectra. 

This article is arranged as follows. First we give some 


background on /3-decay, especially the FF decays, and 
the nuclear structure theories we used in our calculations. 
Then we present the calculated results with comparisons 
to experiments and the corresponding neutrino spectra. 
Conclusions are given at the end. 

II. THEORY OF /3-DECAY 

The decay rate for /3-decay can be written generally as 

0i: 

A = In 2/ti/2 = Ai. (1) 

i 

With the conventions and numerical constants used in 
0,0 one obtains: 

/ = 8896 A 

/ ^Qi 

C{ui)F{Z,uj)puj{ujQi — w^doj. (2) 

I 

Here uj = Ae/rue is the energy of the emitted electron 
in the units of electron mass, wq is the /3 -decay energy 
in the unit of electron mass, and p = \/oF — \ is the 
momentum of the electron. F(Z, oj) is the Fermi factor 
which takes into account of nuclear charge on the shape 
of the spectra for the emitted electron. 

The nuclear structure dependence on the shape of the 
emitted leptons is contained in C(lo). It has different 
w-dependencies for different kinds of decays that lead to 
the different spectra for emitted electron and neutrino. 
For allowed-decay, C{ijo) is independent of w. For the FF 
decay, the dependence can be written in the form Q: 

C{uj) = Kq + Kiu + K_iIu + K2UJ^. (3) 

For FF decays one has three different types of transitions 
associated with the change of spins, AJ^ — 0“, 1“,2“, 
they have different matrix elements and to dependencies: 

C^^=°{oj) = Ko + K_iloj 
C^J=\uj) = Ko + Kiuj + K.i/uj + K2 Uj' 

= Ko + Kioj + K20 j' 


( 4 ) 
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The detailed expressions for the K’s can be obtained from 
0,0- For 0“, there are three matrix elements MJ, Mg' 
and Mq, for 1“ one has five matrix elements involving 
u, u' , X, x' and y, and for 2“ just one matrix element 
0 is involved. The expressions for these matrix elements 
are given in Q. In Q only the Mg, u, x and z terms 
were used for FF branches. Our additional terms result 
in some differences between our results and those of 
To get the electron or neutrino spectra, we take deriva¬ 
tives over the respective energies: 

= C{uj)F{Z,uj)p{uJo - uj)"^ 

auj diu 


dUJiy dOJjy 

= C{ujQ - uJu)F{Z, Wo - Wiy)a;^-\/(wo - - 1 

The spectra for FF decays are different from that of al¬ 
lowed GT, and their shape depends on the decay modes 
(J'"’). To obtain the spectra we need to know some de¬ 
tailed structure information for the /3-active nuclei. 

For the nuclear structure calculations, the configura¬ 
tion interaction model or shell model (SM) provides an 
exact solution within a model space for a restricted set of 
valence orbitals. Realistic shell-model Hamiltonians can 
be derived from renormalized interactions based on the 
nucleon-nucleon interaction with some empirical single¬ 
particle energies and modifications to reproduce exper¬ 
imental binding energies and excitation energies. How¬ 
ever, as the number of valence nucleons increase, the di¬ 
mensions of the configurations increase drastically mak¬ 
ing the calculations impossible. Starting with a closed 
shell of ^^^Sn one add nucleons in the “jj56” model 
space that consists of the five ( 137 / 2 , 2(35/2, 2(33/2, 3si/2, 
lhii/ 2 ) orbitals for protons and the six ( 1 / 19 / 2 , 2 / 7 / 2 , 
2/5/2, 3p3/2, 3pi/2, 1*13/2) orbitals for neutrons. We can 
consider up to four neutrons and four protons in this 
j/56 model space. The SM can be applied to the decays 
of nuclei with both even and odd numbers of protons or 
neutrons. 

To obtain results over a wider region of the nuclear 
chart, one needs to use various approximations. One of 
these is the Quasi-particle Random Phase Approxima¬ 
tions (QRPA) which assumes the excited states of the nu¬ 
clei are small harmonic oscillations beyond the Hartree- 
Fock-Boglyubov (HFB) or BCS ground states. Only two 
quasi-particle excitations are considered in this approxi¬ 
mation. By changing one neutron to one proton or vice 
versa, we obtain the spectra for odd-odd nuclei, this is 
the so-called pn-QRPA method 0 which is usually used 
for charge exchange reactions as well as /3-decay. The 
QRPA method can only be applied to even-even nuclei. 

III. RESULTS AND DISCUSSION 

For the SM calculations we used the NuShellX@MSU 
code [ 1 ^. The Hamiltonian for the j/56 model space 


is taken from m- For this model space, the spin-orbit 
partner of h and i levels are not included. As a result, 
a larger than average quenching is needed for the cal¬ 
culated Gamow-Teller matrix elements. The truncation 
to /j56 will also require renormalization of the various 
FF operators. The goal is to reproduce the experimental 
logit values so that we will have realistic results for the 
neutrino spectra. 

For the QRPA calculations we use a Hamiltonian sim¬ 
ilar to that used in We start with the realistic 
G-matrix for GD-Bonn interaction, and then introduce 
two renormalization parameters, gph and gpp, for the 
particle-hole and particle-particle channels, respectively. 
The fitting strategy for them, as well as for the quench¬ 
ing of factors for both GT and FF decay is explained 
in 0. In this work we generally follow the previous 
work where we used gA(v) = 0-5 9ao(vo) foi' s-H Wpes of 
transitions, but slightly change some of the parameters 
to better reproduce the logit values in the Xe region: 
3^(1, 2“) = 0.43AO, and 3y(0“) = 0.63^0- The same 
quenching values are used for the SM calculations. 

In TableUwe present the comparisons of experimental 
results with the SM and QRPA methods for the two even- 
even nuclei. For ^^®Te, where the experimental data and 
both calculations are possible, we see good agreement 
among them. A one-to-one correspondence of most decay 
branches can be found between the SM calculations and 
the experimental results, the difference of the logit values 
are within 0.2 which means a factor of 1.5 in the transi¬ 
tion rates. The QRPA calculations agree with the shell 
model with differences for logit values around 0.1 — 0.2. 
Another even-even nucleus which has been measured is 
^^°Xe. However, it is beyond the reach of our current SM 
computational capacity, so only QRPA results are shown. 
One finds that for this nucleus, the QRPA calculations 
are in good agreement with the measurement. 

As we have stated above, different decay channels may 
have different shapes due to different dependencies over 
energy w, so we need to investigate the effects of these 
decay channels on the neutrino spectra shape. For the 
odd-odd or odd-A nuclei there is usually mixing between 
different decay channels as \Ji — Jf \ < AJ < Ji + Jf, 
but for even-even nuclei, because the ground states of 
the parent nuclei has always Ji = 0, AJ is unique for 
specific final state of daughter nuclei, there will be no 
mixing among different channels and it is easy to isolate 
different shape changes in different decay channels. 

In figlU we compare the neutrino spectra shape 
changes relative to the allowed shape for different chan¬ 
nels with different methods for two even-even nuclei 
(i36Te and ^^°Xe). For each nucleus we show the 0“, 
1“ and 2~ decay branches. The SM and QRPA methods 
agree well with each other. For 0“ decays, the change 
of the spectra is small and it is a good approximation 
to treat the 0“ decay as allowed decay. For 1“ decay 
the change is large with the peak of the neutrino spectra 
shifted downwards. This means that more neutrinos have 
less energy than expected from the previous simulation^ 
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TABLE I: List of excitation energies and spin-parities of the final states and the corresponding logft values from the experiments, 
the shell-model (SM) and pn-QRPA(QRPA) calculations for different nuclei, we are choosing here only important low-lying FF 
branches. The measured half-lives from [Q are presented here. The excitation energies are in the unit of MeV. For QRPA 
calculations compared with [^, we have minor changes on the quenching (explained in the text) to make it much closer to the 
experimental results in this region for a better comparison. 



J” t(s) 

Jf 

Exp. [12] 
Eex 

log/t 

ShM 
Jf Eex 

log/t 

JJ 

QRPA 

Eex 

log/t 


0+ 17.63 

(1-) 

0 

>6.7 

1- 0 

6.85 

Q- 

0 

6.37 

136 Te 


(0-,l,2-) 

0.222 

7.23 

2” 0.095 

7.37 

1" 

0.171 

6.95 



(o-,i) 

0.334 

6.27 

0" 0.133 

6.41 

2“ 

0.194 

7.89 



(o-,i) 

0.631 

6.28 

1" 0.426 

6.26 

2“ 

0.541 

6.99 



(0-,l,2-) 

0.738 

7.57 

2" 0.507 

6.71 

1” 

0.747 

6.13 


0+ 13.6 

i-,o- 

0.080 

6.14 



0' 

0 

6.15 

i"°Xe 


(0,1-) 

0.515 

6.82 



1^ 

0.127 

6.77 



0(-),i(-) 

0.653 

5.98 



2" 

0.365 

7.01 



(1,2-) 

0.800 

Ri7.1 



1^ 

0.586 

6.05 



l(-) 

0.966 

6.77 



1^ 

1.353 

6.75 



(MeV) (MeV) E^ (MeV) 


FIG. 1: (Color online) Neutrino spectra of low-lying FF decay branches for ^®®Te (left) and (right) from SM (dashed 

lines in la) and lb)) and QRPA (thick colors) calculations. The meanings of different line-styles are illustrated in the graph. 
For the Full FF means this calculations with matrix elements and phase space stated above in the text, meanwhile the 

“Simp. FF” means the simplified matrix element used in [^. 


using the allowed type of phase space. For 2“ decay the 
behavior of the change to the shape is a bit different from 
that of 1“ as seen from figH] where the shape of the neu¬ 
trino spectra for this decay branch is broadened. 

We also make a comparison of the full microscopic 
calculations to the approximations made in Q where 4 
out of 9 matrix elements are used (affecting the 0“ and 
1“ decays). For the 0“ decay, the approximation used 
in Q gives a result that is opposite to the full micro¬ 


scopic calculations, slightly shifting the neutrino spectra 
to lower energy. For 1~ decay, the approximation com¬ 
pletely changes the behavior of the neutrino spectra. Due 
to the over simplified forms in , the behavior of an over¬ 
all shift of spectra to low energies disappears now. This 
comes from the fact that for simplified 1“ decay in table 
I of @ one of its matrix elements ([E,r]^“ or u in this 
work) has the same form as that for 2“ decay ([S,r]^“ 
or z in this work). 
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TABLE II: The percentage of the numbers of neutrinos which 
of the actual decay compared with the allowed shapes used in 
the simulation for single decay branches of ^®®Te and ^^°Xe, 
denoted by S defined in text. The superscripts here are Q 
for QRPA and S for shell model, the subscript “simp” means 
that we used the simplified FF matrix-elements used in [^. 





0“ 

ir 

2r 

0.0 1.002 0.995 0.133 1.001 
0.171 0.899 0.929 0.0 0.902 
0.747 0.938 0.971 0.426 0.933 
0.194 0.968 0.065 0.970 

0.541 0.968 0.507 0.982 

0.0 1.003 0.990 
0.127 0.875 0.949 
0.586 0.919 0.981 
0.060 0.971 

0.365 0.976 


There is similar behavior between ^^®Te and We 

would also expect the same behaviors of these FF decay 
channels in odd-mass or odd-odd nuclei since they have 
the same transition operators as the even-even nuclei. 
From the above results, we conclude that the inclusion 
of FF decays could eliminate the “reactor anti-neutrino 
anomaly” if there are enough beta branches containing 
1,2“ transitions with suitable end-point energies, espe¬ 
cially 1“. However, if we examine the nuclear chart for 
the decay branching ratios, we find that 1,2“ are usually 
accompanied with 0“ decays which usually have a much 
smaller logit values (a stronger transition probability). 
This would reduce the overall changes to the spectra. 

To quantify the change in the neutrino spectrum due to 
the change of phase space, we integrate over the spectra 
with the two phase spaces as follows, 

_ 1 — uppiE < Et) 

1 — ncriE < Et) 

dN 

ni{E<Et) = / —{E,)dE, (6) 

Jo dEt, 

with dN/dEi,{E^)dEi, = 1. Et is the energy 

needed to trigger the interaction Df, + p ^ e~ + n, and 
Ee,nd is the maximum energy of emitted neutrinos. The 
reduction in the number of low-energy neutrinos is given 
by A = 1 — i5. The change depends on the end point en¬ 
ergy Eend, which can be expressed as Qp — me — Eex- So 
we need precise excitation energies for the determination 
of neutrino spectra. This result can then be compared 
with the value of the reactor neutrino anomaly to see if 
the lack of FF phase space factor in the simulation can 
explain the missing neutrinos. The results for single de¬ 
cay branches are listed in Tabldlll A comparison between 
QRPA and shell model shows similarities for the ratio 5, 
this agrees with FiglT] For the detailed values; the change 
A of the 0“ decay is negligible, for 1“, A goes up to ten 
percent, and for 2“, A is only 2-3 percent. 

To obtain quantitative results on the dependence of the 
detailed changes on the end point energies of the decay 
branches, we vary the Q values in the calculations for 
the two nuclei ^^®Te and ^"^°Xe. The results are plotted 
in fig 12] where one observes that to a large extent this 
relation is nucleus independent. A for the 0“ branches 



FIG. 2: (Color online) Dependence of the changes for the 
percentage of neutrino number 5 (defined in text) on the end¬ 
point energies for several FF decay branches of ^®®Te(bold 
lines for QRPA calculations and dashed lines for SM calcu¬ 
lations) and ^^°Xe(dashed-dot lines for QRPA calculation). 
Here we varying the Q values of the two nuclei above to see 
how the changes are related to the end-point energies of the 
decay branches. 

are near zero except below end-point energies of 3 MeV. 
For small end-point energies, A is large due to the shape 
changes at the spectra tail; but these are not important 
since contributions of these branches to the total spectra 
are small, see Fig.3 of For 2“ decay the dependence 
of 5 on the end-point energies are independent of logit 
values since it has only one component. For end-point 
energies from 4-6 MeV, A is around 3 — 4%. 

However, for 1“ decays 5 depends on both Eend and 
logit. To see this we also plot the 1^ decay branches for 
the two nuclei FigI21 Compared to 1)" the FF decays to 
the I 2 states have smaller logit values (TablelU (i.e. they 
are stonger) and have smaller A values (Figl2|). The rea¬ 
son of this comes from the fact that the transition rates 
of 1“ are determined by five different components. They 
are combined to give the final decay rates, and their dif¬ 
ferent combinations have different energy dependencies. 
At Eend ^4 — 6 MeV, A is 5 — 15%. It was estimated 
in @ that 30% of the decay branches of the fission prod¬ 
ucts are FF. Thus, in the most extreme case where the 
FF is dominated by A = 1“ the change of the neutrino 
spectrum could be as large as A = 4.5%. 


IV. CONCLUSION 

In this work, explicit analysis of /3-decay neutrino spec¬ 
tra with inclusion of the first forbidden part has been 
performed. One finds that use of the allowed decay 
phase space factor results in a correction of up to about 
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A = 4.5% due to AJ’^ = 1“ FF transitions. An av¬ 
erage over all types of FF transitions, end-point ener¬ 
gies and logit values would result in a smaller value of 
A = 1 — 2%. The finite size effects and the weak mag¬ 
netism corrections obtained in for the allowed (GT) 
decays are estimated to be A = 2 — 3%. If the average 
branching ratios for all types of FF is estimated, they 
can be combined with our results to obtain an improved 
correction for the shape of the neutrino spectra. 
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